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Abstract 

This study deals with continuous Hmits of interacting one-dimensional diffusive systems, arising from stochas- 
tic distortions of discrete curves with various kinds of coding representations. These systems are essentially 
of a reaction-diffusion nature. In the non-reversible case, the invariant measure has generally a non Gibbs 
form. The corresponding steady-state regime is analyzed in detail with the help of a tagged particle and 
a state-graph cycle expansion of the probability currents. As a consequence, the constants appearing in 
Lotka-Volterraequations — which describe the fluid limits of stationary states — can be traced back directly 
at the discrete level to tagged particles cycles coefficients. Current fluctuations are also studied and the 
Lagrangian is obtained by an iterative scheme. The related Hamilton- Jacobi equation, which leads to the 
large deviation functional, is analyzed and solved in the reversible case for the sake of checking. 
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1 Introduction 

Interplay between discrete and continuous description is a recurrent question in statistical physics, which 
in some cases can be addressed quite rigorously via probabilistic methods. In the context of reaction- 
diffusion systems this amounts to studying fluid or hydrodynamic limits, and number of approaches have 
been proposed, inparticular in the framework of exclusion processes, see [21], [2] [SO], [22] and references 
therein. As far as the above limits are at stake, all these methods have in common to be limited to systems 
having stationary states given in closed product form, or at least to systems for which the invariant measure 
for flnite N is explicitly known. For instance, ASEP with open boundary are described in terms of matrix 
product forms (really a sort of non-commutative product form) , and the continuous limits can be understood 
by means of Brownian bridges [8]. We propose to tackle these problems from a different view-point. The 
initial objects are discrete sample paths enduring stochastic deformations, and our primary concern is to 
understand the nature of the limit curves, when N goes to infinity: how do they evolve in time, and which 
limiting process do they represent as t goes to infinity: in other words, what are the equilibrium curves? 
Following [14] and [15], we give here some partial answers to these questions. 

In [H] a specific model was considered, namely paths on the square lattice, and we could reformulate the 
problem in terms of coupled exclusion processes, to understand the thermodynamic equilibrium and a phase 
transition point above which curves reach a deterministic profile, solution of a nonlinear dynamical system 
which was solved explicitly by means of elliptic functions. Two extensions of this system were introduced in 

m ■■ 

• one which comprises multi-type exclusion particle systems encountered in another context (see e.g. 
|12[ 113)). including the ABC model for which similar features occur [6J; 

• a tri-coupled exclusion process to represent the stochastic dynamics of curves in the three-dimensional 
space. 

With this extended formulation, we provided a set of general conditions for reversibility, by analyzing cycles 
in the state space and the corresponding invariant measure. 

This paper focuses on non-Gibbs states and transient regimes. In another work in progress [H], we analyze 
the asymmetric simple exclusion process (asep) on a torus. Under suitable initial conditions, the usual 
sequence of empirical measures converges in probability to a deterministic measure, which is the unique 
weak solution of a Cauchy problem. The method presents some new features, and relies on the analysis of a 
family of parabolic differential operators, involving variational calculus. This approach let hope for a pretty 
large level of generalization, and we are working over its general conditions of validity. 

Sections [3] and [4] are devoted to the stationary regime, for which, from [I4j and the limit curves are 
known to satisfy a differential system of Lotka-Volterra type which is the essence of the fiuid limits in our 
context. Section [3] solves the steady state regime in the reversible case. A geometric interpretation of the free 
energy is provided (involving the algebraic area enclosed by the curve) , as well as an urn model description 
for the underlying dynamical system, leading precisely to a Lotka-Volterra system. 

Non-Gibbs states are considered in section |4l In [T^ , necessary and sufiicient conditions for reversibility 
where given, by identification of a family of independent cycles in the state graph, for which Kolmogorov's 
criteria have to be fulfilled. We pursue this analysis by showing that irreversibility occurs as a result of 
particle currents attached to these cycles. A connection between recursion properties (originating matrix 
solutions) and particle cycles in the state-graph is found, with the introduction of loop currents, on the 
analogy with electric circuits. These recursions at discrete level connect together invariant measures of 
systems of size N (the number of sites) and of size — 1, and they involve coefficients which are given a 
concrete meaning. Indeed, by means of a functional approach, we map explicitly these structure coefficients 
onto special constants which intervene in the Lotka-Volterra systems describing the fiuid limit, as A^ — *■ oo. 

In the last section [Sj we observe that local equilibrium takes place at a rapid time-scale, compared to the 
diffusion time which is the natural scale of the system. We extend the iterative scheme procedure initiated 
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2 Model definition 



in [14] and developed in [15], which originally concerned only the steady-state regime. In fact, this scheme 
allow us to express in transient regime particle-currents in terms of deterministic particle densities: this is 
a mere consequence of a law of large numbers. At least when the diffusion scale is identical for all particle 
species, local correlations are found to be absent at the hydrodynamical scale. Finally, in the spirit of the 
study made in [3], we obtain the Lagrangian describing the fluctuations of currents, and we analyze the 
related Hamilton- Jacobi equations. 



2 Model definition 



2.1 A stochastic clock model 

The system consists of an oriented path embedded in a bidimensional manifold, with N steps of equal size, 
each one being chosen among a discrete set of n possible orientations, drawn from the set of angles with some 
given origin {^^, fc = 0, . . . , n — 1}. The stochastic dynamics in force consists in displacing one single point 
at a time without breaking the path, while keeping all links within the set of admissible orientations. In 
this operation, two links are simultaneously displaced. This constrains quite strongly the possible dynamical 
rules, which are given in terms of reactions between consecutive links. 



For any n, we can define 



X'^X^ ^ X^X'', k n], k =^ I, (2.1) 



which in the sequel will be sometimes referred to as a local exchange process. It is necessary to discriminate 
between n odd and n even. Indeed, for n = 2p, there is another set of possible stochastic rules: 



^fc^i ^ ^;^fc^ = l^k+p, 

^Ik 

^k^k+p ^ ^k+i^k+p+i^ fc = l,...,n. 



(2.2) 



The distinction is simply due to the presence, for even rt, of folds (two consecutive links with opposite 
directions), which may undergo different transition rules, leading to a richer dynamics. The parameters 
{Afei} represent the exchange rates between two consecutive links, while the 7fc's and 5k s, correspond to the 
rotation of a fold to the right or to the left. 



2.2 Examples 

1 ) The simple exclusion process 

The first elementary and most studied example is the simple exclusion process, which after mapping particles 
onto links corresponds to a one-dimensional fiuctuating interface. In that case, we simply have a binary 
alphabet. Letting = r and X^ = f, the reactions rewrite 

\- 

TT ^ fr, 

A+ 

where A* is the transition rate for the jump of a particle to the right or to the left. 

2) The triangular lattice and the ABC model 

Here the evolution of the random walk is restricted to the triangular lattice. A link (or step) of the walk 
is either 1, e^^'^^^ or e'**'^/'^, and quite naturally will be said to be of type A, B and C, respectively. This 
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corresponds to the so-called ABC model, since there is a coding by a 3-letter alphabet. The set of transitions 
(or reactions) is given by 



AB 



BA, 



BC 



CB, 



CA ^ AC, 



(2.3) 



where the rates are arbitrary positive numbers. Also we impose periodic boundary conditions on the sample 
paths. This model was first introduced in [12] in the context of particles with exclusion, and, for some cases 
corresponding to reversibility, a Gibbs form has been found in [13] . 

3) A coupled exclusion model in the square lattice 

This model was introduced in [14] to analyze stochastic distortions of a walk in the square lattice. Assuming 
links are counterclockwise oriented, the following transitions can take place. 



AB 




BA, 


BC 


Acb 


CB, 


CD 


Ado 


DC, 


DA 


Kd 


AD, 




Kb 




Abe 




Acd 




Ada 


AC 


&bd 
7a c 


BD, 


BD 


Saa 
Ibd 


CA, 


CA 


Sdb 


DB, 


DB 


Sac 

Idb 


AC. 



We studied a rotation invariant version of this model, namely when 



'a+ 


def > 
— ■'^ab 


= \c = ^cd 


= ^da 


A- 


def , 

— Ma 


= ^cb = ^dc 


= Xad 




def 






= 7ac 


— Ibd — lea 


= Idb- 




def r 
= dac 


= 5bd = Sea = 


= Sdb- 



(2.4) 



Define the mapping {A, B, C, D) {t'',t^) G {0, 1}^ such that 



A~ 


^(0,0), 


B - 


-(1,0), 


C - 


-(1,1), 


D - 


-(0,1). 



(2.5) 



The dynamics can be formulated in terms of coupled exclusion processes. The evolution of the sample path 
is represented by a Markov process with state space the set of 2iV-tuples of binary random variables {rf } 
and {tj^}, i — 1, . . . ,N, taking the value 1 if a particle is present and otherwise. The jump rates to the 
right (+) or to the left (— ) are then given by 



A? («) = rf ff+iA± + rf rf+iAT + ff rf+i7± + rf ff+i7^. 



(2.6) 



Notably, one sees the jump rates of a given sequence are locally conditionally defined by the complementary 
sequence. 



3 Stationary regime for reversible systems 



In this section, we quote the main characteristics of the steady state distribution when the processes at stake 
are reversible. 
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3 Stationary regime for reversible systems 



3.1 The general form of the invariant measure 

Up to a slight abuse in the notation, we let e {0, 1} denote the binary random variable representing the 
occupation of site i by a letter of type k. The state of the system is represented by the array rj {X^^ , i — 
1, . . . ,N;k — 1, . . . ,n} oi size N x n. The invariant measure of the Markov process of interest is given by 

TT,, = iexp[-7^(r,)], (3.1) 

where 

^(^)-iEE<'^'^'' (3-2) 

i<j k.l 

with a^^^ and a^^^ two A^-dependent coefEcients related by 

log ^, (3.3) 

provided that some balance conditions hold (see e.g. j2T|)- For example, in the clock model l|2.ip . these 
conditions take the simple form 

EKT'-<')A^. = 0, (3.4) 

and they follow indeed directly from Kolmogorov's criteria (applied to a particle crossing the system), which 
is tantamount to detailed balance equations. 



3.1.1 An example in the square lattice 

To show a concrete exploitation of the form p.ip . we consider the square-lattice model introduced in [14j. 
It does illustrate the rules l|2.2p . Instead of handling the problem directly with the natural set of four letters 
{A, B, C, D}, we found convenient to represent the degrees of freedom by pairs of binary components. In 
the symmetric version of the model defined by l|2.4p . when cycles are absent {Na = Nb = 1/2 and 7^ = 7^), 
we could derive the invariant measure 

7r„ = iexp[/?^(rff^-rff;)], (3.5) 

with 1] = {(rf , Tj^), i = 1 . . . N} with /3 = log Let us see how this relates to the original formulation of 
the model in terms of the four letters A, B, C and D. 

Proposition 3.1. Under the reversibility conditions imposed on the transitions rates {Xki,j^, 5^ ,k = 1 . . . 4, Z = 
1 . . .4}, the measure given by Il3.1]) and ^3.2]) reduces to 

n„ = ^expl!^^B,Aj - A,Bj + A,Dj - D,Aj 

+ C^Bj - B,Cj + D,Cj - Ci^.j, (3.6) 

and is equivalent to ( tg.J)) . 

The proof is not difficult, starting from (|3.5p . It can also be achieved by a direct argument, i.e. without 
using l|3.5p . from theorem 3.2 of [15]. 



3.2 Free energy 
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3.2 Free energy 

We consider again the ABC model as a typical example, and the extension to other models will be straight- 
forward. Assume conditions (|3.4p hold, so that the invariant measure is given by 

N 



1 1 ^ 

= - exp [- aZ'AB, + a<:'S.C, + a^^'C,^,] , 



i<j 

where the constants a^^', a^™' and a^^* take the values 

'^fjo ^cb ^ac 

while a^"' a^^' and ajj"' are set to zero, to be consistent with (|3.3p . The constraints p.4p now become 

Following [6j, we want to write a large deviation functional corresponding to the above Gibbs measure when 
N ^ oo. Set X = jj, J — exp(2i7r/3), and let Z{x) denote the complex number given by 

l^N] .A n n 

t2 



where we have introduced the parameters 

a = lim a^,"', /3 = lim a[Z\ 7 = 1™ a^fc • 

N^co N^oo N—^oo 

The sequence 77 = {{Ai,Bi, d), i — 1 . . . N} is thus represented by a discrete path F in the complex plane, 
made of oriented links having only three possible directions 

{61 = 0, 61 = 27r/3, 6* = 47r/3}, 

depending on whether a particle A, B or C is present. The length of a Hnk corresponding to A, B, or C is, 
respectively, 1/iNa), 1/{N(3) or l/(iV7). 

The equation of F is given by a function Z x Z{x),x G C. Note that condition l|3.7p ensures F is 
closed, that is 

Z(l) = — 1— (1 + J + J2) = 0. 

a + p + j 



The area A enclosed by F is given by 



A'^ — (p [zdz — zdz) , 



and, for large iV, this coincide with 

Vi^AuBk Ck\ , BifCk Ak\ , CifAk Bk^ 



As a result. 



2V3 (q! + /3 + 7) 

The large deviation probability is easily obtained from the law of large numbers. It is given by 

PN{Pa,Pb,Pc) = ^ CXp{-NJ='{pa, Pb, Pc)) , (3.9) 
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3 Stationary regime for reversible systems 



with the free energy 



af3j 
2Vl 



^{Pa, Pb, Pc) - S{pa, Pb, Pc), 



(3.10) 



where 



4if /"^j Pa{x) ( Pb{y) Pc{y)\ , Pb{x)(pc{y) Pa{y)\ , Pc{x)/pa{y) Pb{y) 



A{pa, Pb, Pc) = [ dx [ dy ^"^^^ ( 
Ja Jx o; \ 



P 7 



7 V a (i 



namely the way of 



and where the entropy term comes from a multinomial combinatorial factor — — , 
arranging a box of n=[N dx] sites, with 3 species of identical particles having respective populations ni = 
Pi{x)Ndx, i e {a, 5, c}. Stirling's formula for large N yields 



S{pa,Pb,Pc 



dx[pa{x) \ogpaix) + pb{x) \og pb{x) + pc{x) Iogpc(a;)]- 



Stable and metastable deterministic profiles correspond to local minima of the free-energy. According to 
(|3.10p , an optimal profile is a compromise between a maximal entropy and a minimum of the enclosed alge- 
braic area. Curves of maximal entropy are typically Brownian, and they have an area which scales like 1/iV; 
on the other hand, the opposite extreme configuration consisting of an equilateral triangle with negative 
orientation achieves the minimum algebraic area, but belongs to a class of profiles for which the entropy con- 
tribution is equal to zero (since plogp vanishes both for p = and p — 1). Depending on the ratio af3"i/2^ 
of the two contributions, we obtain either Brownian (the degenerate point of the deterministic equations, 
see below) or deterministic profiles, both regimes being separated by a second order phase transition. 

3.3 Lotka Volterra systems 

Under the scaling earlier defined, letting TV oo, we show on two examples that the limiting invariant 
measure is the solution of a non-linerar differential system of Lotka- Volterra type. 



3.3.1 Urn model 



Consider three species, denoted by {A,B,C}, and let N^'^\t), Nl^'^\t) and N^'^\t) be the corresponding 
time-dependent populations. The system is closed, Na + Ni, + Nc — N. At random times taken as exponential 
events, individuals do meet and population transfer take place at rates a, /?, 7, associated with the reactions 

AB^BB, 

7 



BC^ 

CA- 



>CC, 
■AA. 



This zero-range process is an urn- type model of Ehrenfest Class, as defined in ^8j, where indivivuals, rather 
than urns, are chosen at random. When N increases to infinity, we rather consider concentrations instead 
of integer numbers: 

for i — a,b,c. After a proper scaling limit, the dynamics of the model is described by the following Lotka- 
Volterra system 

-Q^ = Pa[f3pc ~ IPb), 

dpb I X 
= PbKlPa - aPc), 



dpc 
dx 



Pc{apb - (3 Pa), 



3.3 Lotka Volterra systems 
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which, after replacing xhy t and densities by concentrations, is nothing else but the differential system giving 
the invariant measure of the {A,B,C) model, in the fluid limit at thermodynamical equilibrium [6]. 



3.3.2 The square lattice model 

From l|3.6p . one can write down the large deviation functional J^{pA, Pb,Pc, Pd), [as in l|3.9p ]. together with 
the conditions ensuring an optimal profile. This leads again to a differential system of Lotka- Volterra class 

dpA , s dpB , . 

^nPAKPB ~ PD), ^r\PB{pc~ PA)-, 

dpc I ^ dpD , . „ . 

-g^ ^ VPc[PD - Pb), ^t]pd[pa- PC), (3.11) 

in which the last equation follows merely by summing up the three other ones. This system is structurally 
different from the one obtained in [14], which involved only two independent profiles {pa,Pb) corresponding 
to deterministic densities for the particles and while in the present case there are three {pA, Pb,Pc for 
example) . 

It is interesting to notice that, in both models, expHcit level surfaces exist. Indeed, the above system satisfies 
PaPbPcPd = cte, in addition to constraint Pa + Pb + Pc + Pd = 1- On the other hand, Pq(1 — Pa)pb{^ — Pb) 
is the level surface of the former system encountered in [14]. This can be explained by reversing the mapping 
((3131) . so that 

At — , Oi — Tj^ , 

a = , A = ff T^ (3.12) 

This indicates that the set of 4-tuples {"rf , ff , rf , fj^} constitutes the elementary blocks of the system, and 
that letters Ai, Bi,Ci, Di are composite variables encoding correlations of these building blocks. Therefore, 
in the continuous limit, we are left with two different descriptions of the same system, related in a non trivial 
manner. We propose now to explore more carefully this connection. In particular, while the linear mapping 

+ (3.13) 

still holds in the continuous limit, as a relation between expected values 

U = PB+PC, (3^^^^ 
[Pb ^ PC + PD, 

the non-Hnear equations (|3.12p are instead expected to bring a different form, since they involve correlations. 
Proposition 3.2. The differential system given by 

is related to 113. through the invertible functional mapping given by 

ipA^PaPb + K, pB=PaPb-K, 
\PC = PaPb + K, pD = PaPb - K, 




where K is a constant to be determined. 
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4 Non-Gibbs steady state regime 



Proof. First, let {pb, Pc, Pd} be the set of independent variables in l|3.1ip . and express them in terms of 
the new triple {pa, Pb, Pc} given by (|3.14p . This gives 



djPa - Pc) 
dx 

d{pb - Pc) 



dx 



dp, 



'c 



dx 



fliPa - Pc){Pa + Pb- 1), 
ViPb - Pc)(l - Pa+ Pb), 
VPciPb ~ Pa)- 



(3.17) 



Combining these equations yields 

(dpa 

dx 
dpb 
dx 

which in turn allows to express pc as 

PC 



VPaiPa + Pb-l)+ Wc(l - 2pa), 
VPb{i - Pa- Pb) + VPcC^Pb - 1), 



(3.18) 



dpb 

Pa~ Pb\ dx 



Pb- 



dpa 

dx 



Instantiating this last value of pc in (|3.18p and in (|3.17p . we obtain (|3.15p . after immediate recombination, 
together with the relation 

dpc _ djpaPb) 
dx dx 

This last equation has its counterpart for pA, Pb and po- after integration, we are left with four constants, 
which reduce to the one given in l|3.16p only when compatibility with p.l4p is imposed. 



4 Non-Gibbs steady state regime 

We call non-Gibbs steady state regime, a regime for which the invariant measure is not described by means 
of a potential. This occurs when reversibility is broken, that is when there exists at least one cycle in the 
state graph for which the Kolmogorov criteria fails. A complete set of detailed balanced equations cannot 
be written in such a case, there exist at least two states rj and 77', connected by a single particle jump, with 
rate A^^' , A^',, such that 

Arjrj'TT,, - A,,/^7r,,/ = (/) 7^ 0, (4.1) 

if TT^ denotes the invariant measure. It is the second member of this equation we wish to study in this 
section. In the sequel we note S the state space, G the corresponding state graph, by assigning oriented 
edges between pair of nodes (a,/3) € S^, when the rate Xap is non-zero, C will denote a cycle in G and we 
denote T the set of spanning trees on G- 



4.1 The tagged particle cycle 

Cycles in the state graph for the n odd model are important in the analysis of reversibility, and they are 
the ones for which at least one particle performs a complete round-trip. For example if a given particle 
makes — 1 successive jumps to the right, because of the circular geometry, the initial and final states are 
identical, up to a 1-step global shift to the left of the particles. As long as this particle is the only one in 
movement, the permutation order of the remaining other — 1 particles is kept frozen. The corresponding 
subsequence ry<""^' will in the sequel denote these specific cycles. Let us examine this one particle model, 
by tagging a specific particle which is given a new label Y, and by following its motion conditionally on 



4.1 The tagged particle cycle 
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ij* — {X^, i — 1 . . . N, k E {1 . . . n}}, the complementary frozen set of particles. This is equivalent to consider 
Y moving in the inter-sites {i + 1/2, i = 0. . . A'^ — 1} of the iV — 1 frozen particles. The question is then 
to analyze the steady-state regime of a particle moving around a circular lattice in a random environment. 
To all allowed transitions which are jumps of Y between sites » — ^ and « + 5, we let correspond the set of 
conditional transition rates given by 



n 



k=l 



Violation of condition l|3.4p corresponds to have, 



N-1 



i=0 i=0 

This coefficient attached to the cycle rj* is the determinant of the set of flux equations 

A+(zK_i -A-(zK+i =(/.(7?*), z-0,...,iV-l, 
giving the invariant measure 7rj_|_i which reads, 

N ^ n,N 

'^^+1 = 7 E E E log ^yr- + E 



(4.2) 



(4.3) 



m=l l + l<j<i 



i<j<l 



0,. 



where Z is a normalization constant. A diagramatic representation of each term in the summation (over 
I) is given in Fig. 14. li b. Each term is in fact a spanning tree on the reduced tagged-particle state-graph, 
weighted by the transitions rates and rooted at the considerd point {i + h for tt. 



The constant Z is 

therefore the sum of all spanning-trees on the reduced tagged-particle state-graph. The probability current 
between site i — h and site i + k reads 



exp(^ N^logXyr, 



exp(^ 7V„logA 



with Nm the number of particles of type to, a quantity independent of z. This shows that (j){ri*) is a quantity 
attached to the cycle rj* , which will be referred to as cycle current and reads 



0(^*) = |det(7y*). 



(4.4) 



Depending on the sign of det(?7*), the diffusion of particle Y is biased in the right {Aet{ri*) > 0) or in the 
left (det(77*) < 0) direction. Of course the reversible case is recovered when the determinant vanishes, which 
corresponds exactly to Kolmogorov's criterion. 

Case of open systems: example of ASEP 

Consider the well studied asymmetric simple exclusion process ASEP with open boundary conditions, defined 
by a the rate of particle entering to the left side and f3 the rate at which particles exit from the right side. 
The generaHzation to open systems of our definition of the tagged particle cycle (tpc) is depicted in figure 
14. 21 a. We adopt the convention for the cycle orientation that particles move positively to the right and holes 
to the left. Assume we give a tag to one of the particles. Let it perform sucessive jumps until reaching 
the right side; when it leaves the system it is in fact transformed into a hole; We keep the tag attached 
to the hole which performs successive jumps in the opposite direction until it reaches the left side; again it 
transformed back into a particle which in turn performs jumps to the right until the reaching of the initial 
position, to conclude the cycle. 
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4 Non-Gibbs steady state regime 




(b) 



Fig. 4.1: (a): relative motion of the tagged particle, (b): corresponding state space and a spanning 
tree contribution to tts 

o»oo»» o 



o»o o»o « 




o»oo»»» 



o • o»o»« 



oo» oo •• 



»o«oo •• 



Fig. 4.2: Example of a tagged particle cycle in the state graph for asep with 7 particles and open 
boundary. 



4.2 Combinatorial formulas for invariant measure and currents 

We give here a combinatorial way of expressing the stationary measure on a connex finite state space S 
of size J\f, the number of states. We consider a continuous-time irreducible Mark;ov chain, with a set of 
transition rates Xap between states a and /? and define the corresponding state graph G, based on S, by 
assigning oriented edges between pair of nodes (a,/?), when the corresponding rate Aq/3 is non-zero. 

Proposition 4.1. The invariant measure iTa is given by 



(4.5) 



where T is the set of spanning tree over Q, Ta is the set of spanning tree over Q rooted in a, and w(t) the 
weight of a tree t given by 

W{t) = Y\_ ^afj- 

Proof. This follows from reexpressing the solution to the steady-state equation 

where G is the generator, and Gap = "(j^-y ^a-y^Sap + \ap, using the Cramer relation. Indeed since 

Gap = 0, 

/3=1 



4.2 Combinatorial formulas for invariant measure and currents 
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the set of steady-state equations is of rank N — 1 and tt^ can be written as the ratio of two determinants, 
namely the cofactor GaM of Ga^^ and the determinant IGj of the matrix obtained from G by replacing G^ja^ 
by 1 for (3 = \,...,M. G has a structure of an admittance-matrix, as a result, expanding Gap and \G\ 
amounts to sum over spanning trees, 

which leads to formula l|4.5p . ■ 



From this observation we deduce a way to express the probability currents at steady-state, which generalizes 
formula l|4.3p and (|4.4p . First call det(G) a coefficient attached to each cycle G, 

det(G)= n ^-z*- n ^^7. 

(7,<5)ec {'r,s)ec 

generalizing l|4.2p and where the orientation of G is prescribed by the orientation of (a, /3) and the product 
over the set (7, S) e G, is understood according to this orientation. Let Cap the set of cycles in Q containing 
the oriented edge (a,/3). Let 7c a set of subgraph oi G, s.t. when C is glued into a single node ac, Tc 
represents the set of spanning trees rooted in ac- 

Lemma 4.2. The steady state current between states {a, (3) G is given by 

XcpTTa - XpaTTp ^ ^ ^^^^^^^ dct (G) (4.6) 
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Fig. 4.3: (a): state-graph with N = 8 states. Arrows indicates possible transitions, (b): a contribu- 
tion to ttq. (c): a contribution to tt^. (d): a combined contribution to J14. 

Proof. When tTq, is multiplied by Xap, each spanning tree contribution is transformed by drawing an oriented 
edge between a and (3. Since the spanning tree contains by construction of tTq a path going from (3 to a, the 
added edges contributes to the forming of a cycle which contains a and f3. If each oriented edge in this cycle 
have a reversed counterpart, then in A^qTt^ there is to be found a corresponding term with the same edges 
but with reversed orientation in the cycle (see Fig. 14. 3p . In any case, det(G) factors out of an ensemble of 
contributions which consist in drawing trees spanning all the subgraph Q with end-points on G, divided by 
the global normalization constant X^ter ^(^)- This complete the justification of formula l|4.6p . ■ 



Note that J2teTc ^(^) represents the unormalized invariant measure of ac on the reduced graph 

G/C. This indicates that l|4.6p bears recursive properties which could be used for asymptotic limits when 
the size of the system tends to infinity. Let us call G a reversible [resp. non- reversible] cycle if det(G) = 
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4 Non-Gibbs steady state regime 



[resp. det(C) 0]. In the loop expansion of the currents provided by l|4.6p . only non- reversible cycles do 
contribute. For particle system, this distinction is embedded into a topological classification of cycles with 
respect to their corresponding determinant value det(C). 



Connection with the matrix ansatz for ASEP 

For the ASEP model, a simple algorithm has been discovered [9] to obtain the steady-state probabilities of each 
individual state with the help of a matrix ansatz . In this representation, a given sequence i] = 1010 ... 00 is 
represented by a product of matrices D (for 1) and E (for 0), and the corresponding probability measure is 
obtained by taking the trace 

TT,, = ^Tv{WDEDE . . . EE) , 

where W is an additional matrix which takes into account the boundary property. A sufficient condition for 
this to be the invariant measure is that Z?, E, W satisfy 



XiqDE - XoiED 
DW 



D + E 
1, 



(4.7) 



W 



WE = 



iw. 

a 



If Aoi — 0, the process is totally asymmetric (tasep), particles can jump only to the right. Consider the 
system with only 3 sites, which graph is depicted in figure FOl Using these rules we find e.g. that 



1 3 



TTlOO 



.1 

■ a 



-)- 



1 



(4.8) 
(4.9) 



Comparison with the spanning tree expansion is done by counting deletions. A spanning tree is obtained 
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Fig. 4.4: (a):Graph of the state space for a TASEP with three particles and the dual graph corre- 
sponding to the possibles cycles, (b) spanning tree contributions to vrooo- 
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from the complete graph by the deleting of a certain number of edges, and each deletion is accounted for by 
dividing with respect to the corresponding transition rate. The set of spanning trees contributing to ttooo is 
given in figure 14741 b. The brut result (without normalization is): 



/I 1 1\ 



The factor + ^ + /Jl) shows up for each state, and disapears after normalization. Nervertheless, it 
induces in this simple example a factor of 3 in the enumeration of terms, by comparison with the matrix 
ansatz. An underlying symmetry of the state graph is at the origin of this combinatorial factor. Indeed for 
the ASEP system, the steady-state probability current between two sequences -q and -q' separated by a single 
jump between site i and i + 1 reads, 

Aio7r< - AioTT^, = tt;. + 7r;._^^ ^ (4.12) 

as a consequence of l|4.7p . with the subsequence r/* [resp. rj*^i] of rj obtained by deleting bit i [resp. i + 1]. 
We have not been able yet to fill the gap between l|4.6p and (|4.12p . We believe that the combinatorial 
arrangement which occur is due to a hierarchical structure of the state-graph, revealed with the help of the 
tagged particle. The complete analysis of (|4.6p is the subject of another work in progress. Beforehand, in 
the next sections, we simply propose a possible general form for the detailed current equation (|4.ip . which 
leads (see section [474|) to the correct form of the Lotka-Volterra equations describing the fluid Hmits at steady 
state. 



4.3 Cycle currents 

We interpret relation (|4.12p in terms of cycle currents. A transition taking place between two particles of 
different types, say AB — > BA, can be viewed either as a particle A travelling to the right or, conversely, 
as a particle B travelling to the left. In this exchange two joint TPC are involved. In the state-graph, each 
TPC deflnes a face, which we will identify with a subsequence ry*, obtained from rj by removing the tagged 
particle. Accordingly, we attach a set of variables {0(?7*)} G M to each TPC face, while currents between 
states are variables attached to the edges of the graph. Conservation of probability currents at a given node 
is automatically fulfilled, provided that if one write (assuming a transition between site i and i + 1), see 
figure [431) . 

AafcTT^ - A&aTI",,' = 4>a{Vi) ~ 4'b{fli+l) ' (4-13) 

which is tantamount to changing current variables into cycle variables. 

The right-hand side members in l|4.7p and (|4.12p is reminiscent of the second member of (|4.13p . In fact we 
have 

MV*) = TriWrj*) 

With each edge of the state-graph, we associate such an extended detailed balance equation. Then, elimi- 
nating all (^'s from this set of equations leads to the invariant measure equation. Consider the example given 
in figure l4?5l The transition rules are 

AB^BA AC^CA BchcB. 

q 

The various weights corresponding to each sequence and subsequence associated with cycles are given in the 
following table, for g = and g = 1. Note that one should expect tTci = \ and 7rc2 = § from the subgraph of 
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figure 1431 The correction results from the different degeneracy w.r.t circular permutation symmetry (4 for 
Ci and 2 for C2). 
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In this two cases one has the decomposition of the dual variables <f) of relation 14.131 according to 



<j)x{ri*) = C^-K^. with X G {a, 6, c} (4.14) 
with the value of the structure coefficient also given in the table. For example we have 

TTl — -77r2 = CbTTbl ~ CcTTcl- 

This decomposition is however not valid for arbitrary q. A certain number of compatibility constraint have 
to be imposed on the (j)' , because the TPC do not constitute a complete bases of cycles in the state graph. 
When considering the complete system (|4.13p of detailed currents, we have at hand m equations, m being 
the number of edges of the state-graph, and n + vtpc unknowns, where n is the number of nodes and vtpc 
the number of TPC. In matrix form, this reads 

MU = (4.15) 

where 

• M is a m X n matrix; 

• Ha column vector of size n, with the elements the invariant probability measure; 

• $ is a column vector of size m, where each component I is the algebraic contribution of the (two in 
general) TPC having the edge corresponding to I in common. 

To fix the sign conventions, we agree that orientations of cycles are given by the natural orientation of the 
system, i.e. each particle travels positively from left to right. An exception is made for the simple exclusion 
system, since in this case holes travel positively to the left and there is only one type of TPC. 

AioTT,, - AoiTT,,, = <f>{f]*) + 0(?7*+i) for asep. 



AabTT,, - AfcaTT^' = MVi) " MVi+i^ ^^"^ multi-typc systems. 



4.3 Cycle currents 
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(with + 1) the sites involved in the transition). From basic graph theory (see [2]), the quantity giving 
the number of independent cycles in an arbitrary graph G is called the cyclomatic number 



where n,m and p are respectively the number of nodes, edges and components. In our cases, the system is 
irreducible, so p — 1. Since m is the number of equations and n + ^'tp(, the number of unknown, the system 
is over-determined by a quantity 

This over-determination is understood as follows. To each line of the matrix M corresponds a transition 
between two states, so that a given cycle in the state-graph corresponds to some combination of lines of M 
(namely the successive transitions taking part in the cycle), and the resulting sub-matrix is a square matrix 
of size the number of states visited by the cycle. The corresponding determinant vanishes for all trivial 
cycles. Hence the number of independent equations is to — v+vtpc, which is equal to the number of unknown 
minus 1, the remaining degree of freedom being related to the global normalization condition. However, a 
certain number of compatibility conditions have to be imposed on the 0's in order to eliminate safely all 
dependent equations of our system l|4.15p . These conditions are somehow related to the basic recurrence 
scheme which is at the origin of matrix-solutions obtained in the context of ASEP, but also for multi-type 
particle systems [1]. Let us see how the specific form l|4.14p encountered precedingly do combine with these 
compatibility conditions. 



Lemma 4.3. The form 



of the cycle currents fulfills the compati 



/^(N)^(iV-l) 



c. 



(N) (N- 



(4.16) 



imposed by trivial cycles if and only if 
ya,be {1 ...n}. 



Proof. Instead of proving this for an arbitrary trivial cycle, we do it for the one depicted in figure 14. 6[ 
the completion of the general case follows by recurrence, since any trivial cycle can be constructed as a 
combination of cycle of this type. To fix some notation, let 77^, 77^, 77'^ and rj'^ be the states visited by the 



.AB...DC 



.AB...CD. 



.BA...DC . 



.BA...CD. 



Fig. 4.6: Example of a reversible cycle. 



cycle, with i the position of A and j the position of C in rji, so that 
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The sub-system of 14. 151 corresponding to this cycle reads 





— Aba7r^2 ~ 


Mvll- 




(a) 




~ ^dcT^Tj^ — 






(b) 




— ^abT^r)^ — 


Mvll - 




(c) 




— AcdTT^i = 


Mvp] - 




id) 



As already noted, these equations are not independent. Hence the combination Xcd{a)+Xbaib)+Xdc{c)+Xab{d) 
ehminates one equation, but with the resulting constraint on the (/)'s: 

XcdMvl*] ~ XdcMvf+i] + >^dcMvf*] - XcdMvl+i] + 

XbaMvT] - ^-bMvj*+i] + ^abMvT] - XbaMvT+i] = 0- (4-17) 

and rif^i are in correspondence through the transition CD DC at site j, j + 1, as well as rjj* and 
with respect to the transition AB BA at site i, i + 1 .... From the hypothesis of the lemma, (|4.17p 
rewrites 

^(«, /^(«-i,^(«-2, _ ^(«-i)^(«-2)\ +c/("' fci"-^'7r<l?' - Cr^''^'!?') = 0, 

where 77^^** is the sequence obtained from 77^ by suppressing letters at site i and j.The elimination of letters 
in sequences is a commutative process, therefore this last equality holds because of the following identities: 

„1** _ „2** 3** _ 4** 2** „ 3** 4** _ 1** 



The complete study to establishing the range of validity of the recurrence relation (|4.13p altogether with 
(|4.16p is the object of another work in progress. We expect that in general this relation to be valid only 
asymptotically for large N, which could be proved possibly by selecting the dominant terms in the expansion 
631). 



4.4 Fluid limits 

In this section we examine how the microscopic coefficients C^"' , whenever l|4.16p holds, can be transposed 
at macroscopic level and how they are related to important coefficients showing up in the Lotka-Volterre 
equations of the fluid limit. Using the preHminary study ^6\, where a new functional method was introduced 
to handle the hydrodynamic limit of a simple exclusion process, we consider hereafter the n-type case. 



4.4.1 Functional approach 

Let fc = 1 . . . n a set of arbitrary functions in C^[0, 1], G*"' Z/NZ the discrete torus (circle). For 
i e C*', X^{t) is a binary random variable and, at time t, the presence of a particle of type k at site i is 
equivalent to ACf (t) = 1. The exclusion constraint reads 

n 

fc=i 



4.4 Fluid limits 
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The whole trajectory is represented by /^'"'(i) = {Xf'{t),i G G'"',A; — 1 . . . n} which is a Markov process. 
ri<"' will denote its generator and JFj"' — ct(7]<"'(s), s < i) is the associated natural filtration. 

Define the real-valued positive measure 



z(N) rn def 



exp 



where 4> denotes the set {(/)k, k ~ 1 . . . n}. In [16] the convergence of this measure was analyzed for n = 2. 
A functional integral operator was used to characterize limit points of this measure, these were shown to be 
indeed the unique weak solution of a partial differential equation of Cauchy type. 

In what follows, we will be interested in the quantities 

9r{cf>)'^'\og[E{zrm^ 

respectively the moment and cumulant generating function. The idea of using is that the generator, 

when applied to Zl"\ can be expressed as a differential operator with respect to the arbitrary functions (j). 
Indeed, we have 



with 



after having set 



A.(.iV)^-'2A.(A^)e^sinhf^'^'='^^) 



V 2N 

We introduce now the key quantities for hydrodynamic scalings, by assuming an asymptotic expansion of 
the form 

Xki{N)^D{N^ + ^N)+0{l), ykjk^l, 

where aki = —aik are real constants. Here the system is assumed to be equidiffusive, which means there 
exists a constant D such that, for all pairs {k,l), 



lim 



XkiiN) 



D. 



From now on we will omit the argument of Xki{N) and retain the initial notation Xki- The coefficients aki 
express the asymmetry between types k and I. Now one can write 



dt 



(4.18) 



To rearrange the sum in (|4.18p . in order to select dominant terms in the expansion with respect to we 
make use of the exclusion property, which is formally equivalent to 



E' 

fc=i 



d_ 



1 

N' 
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Since we are on the circle i G G'"', Abel's summation formula does not produce any boundary term, so 
that, skipping details, l|4.18p can be rewritten as 



dt 



■ I 1 



I 



r dfr 



dfr 



,2 f(N 



l^k 



It is worth remarking that operators like 



dMjj) 50, (^) 
d'fi 



,2 f(N) 



d 



(4.19) 



d 



d , d 

and 



d 



= ( at) 



and (j^ki^^) — 4>k{jf) produce a scale factor l/N, while 



the above expansion. 



scale as 1/N'^: this explains the selection of dominant terms in 



Let iV — > cxD and assume the convergence of the sequence /q"'. Then, from the tightness of the process, 
together with a zeste of variational and complex variable calculus, as in [16j, we claim [the proof is omitted] 
/i"' also converges, in a good tempered functional space, and its limit ft satisfies the functional integral 
equation 







dx 



d dft 
dx d(j>k{x) 



d'ft 



fe=l ^ ' ' l^k 

Similarly, the cumulant characteristic function is a solution of 

r-l 

= n / 

dt 



d(l)k{x)d(l)i{x) 







d dgt 



k=l 



dx 



dx d(t)k{x) 

dgt dgt 



d^gt 



l^k 



d4>kix) d(l3i{x) d4>kix)d4>i{x) 



(4.20) 



Assume at time the given initial profile pk{x,0) to be twice differentiable with repect to x. Then l|4.20p is 
given by 



9t 



»1 n 

(0) = / dxy2Pkix,t)(l3k{x), 
•^0 fe=l 



where pk{x,t) satisfy the hydrodynamic system of coupled Burger's equations 

k — I, . . . ,n, 



d^Pk ^ 
dx^ dx V 



aikPkPi 



l^k 



with a set of given initial conditions Pfc(a;, 0),k — 1, . . . ,n. 

Remark It is important to note that, without differentiability conditions for the intial profiles Pk{x, 0), one 
can only assert the existence of weak solutions (in the sense of Schwartz's distributions) of Burger's system. 



4.4.2 Functional equation at steady-state 



Theorem 4.4. Consider a particle system of size N, with rules \2.1\ with n types of particles and periodic 
boundary conditions. Assume the detailed current equations holds, for any pair of particle types k and I, 



^kl "V ^Hk "V 



' 'li+l 



kA^l. 



(4.21) 
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Then the limit functional foo [<f>] = lim f^-' [</)] , wh 



.N 



fc=l,i=l 



satisfies the equation 



dx d(j>k{x) 



Ed'^ foo 



l^k 



d(j)k ix)d(j)i (x) 



Cfc/o. 



df. 



d<j)k{x) ' 



(4.22) 



under the fundamental scaling 



Jim log = aki and \fl k, lim J: 



lim — ^ 

N-*oo D 



with 



def \ ^ 



Proof. We use the notation of section 14.41 In order to extract additional information at steady state, we 
refine our preceding variational analysis by defining the functional 



.N 



(«) 

kl 



+ x 



(N) 



92 



r(N) 
I 00 1 



(4.23) 



which corresponds to the second member of equation (|4.18p at steady state, and where it is understood that 

the sets {(f>} {0( ]^), i = 1 • ■ ■ and {d<f>} {^(jf)} « = 1 ■ • ■ N} are taken as independant parameters. 
This functional can be writen in two different manners. Recalling the definitions 

At^kiit) = + 1) - M^) - 0i(« + 1) + M^) = A^fc(i) - AM^), 

(|4.23p may be rewritten in the form 



fe=i.i=i 



(AT) 



5/ 



(AT) 



2 



2 AN) 



d'f. 



- _,(^)a0,(i±i) dM^r)9Uji) 



o{—) 



(4.24) 



On the other hand, combining the sums in (j4.23p yields 



..N 



fc,Z=l,i=l {,,} 



2N 



\(N) (N) UN) (N) 

\l ~ ^Ik ^Tm 



(4.25) 



where 77 is a given configuration, TiX] being the one obtained from 77 by exchanging i and i + 1, and the 
shorthand notation 



k=lA=\ 
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From the assumptions in the statement of the proposition, we can rewrite l|4.25p as 



T<"'({0,9,0}) =iV2 ^ ^e^'^-''+^^^'='«sinh 
k,i=i,i=i {ij} 



2N 



^(N) (N-l) ^(N) (N-1) 



(4.26) 



where rj* is the sequence obtained from rj by removing the site i. We also have 

{';} 

where f^^^''[4>*] means that /^"^' is considered as a function of the n{N — 1) variables {0fc(-^),fc ~ 
1,. . .,n;j = 1, . . . ,N,j ^ i}. Using all these ingredients, expanding (|4.26p in powers of jj and keeping the 
dominant terms, we get 



k,l = lS=l 



(4.27) 



Now, rearranging the summation, using the exclusion property 



d 



1 



comparing (|4.24p and l|4.27p . we finally obtain 



k=lA=l 



■kijj) dM'-^) 2 Va0fc(^)50,(^) dM'-^)dMjj) 



2fW 

oc 



n,N ^^(N) n ^(N) Qi-(N-l) 



As the last equality holds for any dx4>k, letting N ^ oo implies easily (|4.22p . which was to be proved. 



4.4.3 Lotka-Volterra systems and out-of-equilibrium stationary states 

Here we will make the link between the structure coefficients of the current equations (|4.2ip and the fluid 
limit description of stationary states. A solution is sought of the form 

N 



.1 iv 

/oo(0) = expf / y]p^(x)(/<fc(x) 
^•^0 k=l 



which, instantiated into (|4.22p . yields gives the following equations for the p^'s . 

-E^-p^Y.a'^'p^^Ck-vpl^, fc = l...n. 

l^k 

The interpretation of this system is now quite clear : it is exactly a particular stationary solution of the 
system formed by the coupled Burger's equations 



dpk _ d^pk _ d_ 
dt dx^ dx . 

l^k 
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where the functions pk are sought in the class 

Pk{x,t) = p'^{x-vt), 

the variable {x — vt) being taken [modulo 1]. Hence, there is a frame rotating at velocity v, in which is 
periodic. Moreover, in this frame, the stationary currents do not vanish and have constant values 

i^k 

Therefore, while the macroscopic constants {c^, k = 1, . . . , n} are in principle determined from the periodic 
boundary conditions constraints and from the fixed average values of each particle species, they can also be 
directly derived from the microscopic model. 



5 Transient regime and fluctuations 

The goal of this section is twofold : first, establish relationships between currents and particle densities at 
the deterministic level by means of the law of large numbers; secondly, compute the stochastic corrections 
to these relationships for large but finite systems by using central limit theorems and large deviations. 



5.1 Time-scale for local equilibrium 



In keeping with our approach, we discuss the question of local equilibrium by means of the following 
functional 



def 



f = exp 



k,i=i.i=i 



Without entering into cumbersome technical details, let us just notice that the explicit computation of 
L<"'r/"' shows that L<"'r/"' scales like 0{N) instead of 0(1) as L<"'Z<"'. This fact can be interpreted as 
follows. The empirical measure 

n,N 

k,i=i,i=i 

is a convolution of the distribution of interfaces between particle domains with a set of arbitrary functions. To 
any given particle density distribution, drawn from the set of local hydrodynamic densities, there corresponds 
an arrangement of these interfaces which somehow characterizes the local correlations between particles. At 
steady-state, at least in the reversible case, it is easy to show that these correlations vanish. Moreover 
this scaling tells us that correlations vanish at a time-scale faster than the diffusion scale, by a factor of 
N. Therefore, even in transient regime, we expect correlations to be negligible for the family of diffusive 
processes under study. A more formal proof of this fact is postponed to the completion of the functional 
approach initiated in [16] . 



5.2 Hydrodynamical currents 

In our preceding studies, we devised a scheme to obtain a fiuid limit at steady state, first for the reversible 
square-lattice model in [14], and also for the non- reversible ABC model [15]. Here we generalize this procedure 
to transient n-type particle systems, resting upon the hydrodynamic hypothesis, which will be precisely 
stated. The principle of the method is to reverse the relationship between particle and current variables in 
a suitable manner, in order to apply a law of large numbers. 
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5 Transient regime and fluctuations 



5.2.1 DifFusion models 

The system corresponds to rules (|2.ip . For any particle-type k, the rescaled discrete current reads 

Jr (^) + - A,T(^)X,'Vl, z = 1, . . . , iV, (5.1) 

with 



TV 



where arbitrary constants Tk have been introduced (they not modify the value of Jk) to ensure that the A^'s 
never vanish. To be consistent with other scalings, Tk is assumed to scale like N. Our hypothesis is that Jk 
has a limiting distribution, Jk{x), such that, for any integrable complex- valued function a, 



N . . 1 

lim ±ya(A)j(«)(A) ^ / 



a{x)Jk{x)dx. 



(5.2) 



In addition, the system will be said equidiffusive, if there exists a single diffusion constant D, such that, for 
all pair of species (fc, I), 

Xki{N) 



lim 



= D [cquidifFusion] . 



To simplify the notation, consider equation for fc = 1, writing Ja == Ji and replacing Xf by Ai. Then solving 
(|5.ip as a Hnear system yields 

A+(z + l)A.-Jr'(i^) 



Ai+1 = 



This relationship between Ai and Ai+i can be iterated, by means of a 2 x 2 matrix products. Indeed, 
introducing the pair of numbers {m, Vi) such that Ai = the recursion becomes 



Ui+l 
Vi+l 







Ui 




Ui 












yi_ 



where for convenience we divided everything by the common factor y Aa {i)\t{i + 1). Let us define the 
matrices {p being a positive integer) 



^0/ » +P * N def 



G 



■ i+P j_^ def 

■ N ' N' 



S(-) = 



n 



i+P 
3=i 



-- 



Aj(i+1) 



VA+(i+l)A-(i) 
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(explicit references to the species (a) and the size N is omitted here, to Ughten the notations). Because of 
the upper triangular structure of S, we may simply express G as 



N ' N 

To handle this equation in the continuous limit, we need an additional transformation. Define 

. / ^»(') 







Ri 











Ao(») 



together with 



^ TV ' TV 



■ i + p i 



-)R,. 



N ' TV^ 



' iV ' TV^ 

Then G, G° and S verify the same relation. 



(5.3) 



TV ' iV^ 



AT ^N' ' iV '^^^>^ y 



N 



N ' TV 



(5.4) 



but 



S(-)^ 







A^(i+l)A-(i) 

Lo 

Noting that Ai+p+iT a / \a ("^ + P + 1) — ^i+p+i and v4ira/AQ;(i) = Ai, the iteration between i and z +p gives 

Gii(^, ^)A, + G'i2(^, ^) 



A, 



i+P+i 



G9,(^ 



(5.5) 



We can now take advantage of the law of large numbers in equation l|5.4p . First of all, for TV large, and fixing 

we have, 



— i/N and y — _p/TV, letting a = 



-1 



j=i+l,fe=2 fc=2 

from the hydrodynamic hypothesis. To proceed further, we have to distinguish between two situations. 
[The equidifFusion case] 

Recalling that Ta is a free parameter which scales like TV, it is convenient in the equidiffusion case to impose 
the limit 

lim £^ . D. 

N^ao TV 

Then, expanding T,{i/N) with respect to 1/TV yields 



ND 




0(TV- 
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and the limit 



is provided by equation l|5.4p . Hence 



with 



x+y 



Gix + y,x) = Q°{x + y,x) + / du g°{x + y,x + u) :^{x + u) g°{x + u,x), 



G°iy,x) = ex.p\^- du^aakPk{u)j, and 
•^^ k=2 



E{x) 







.7("'(:r)' 







(5.6) 



(5.7) 



still by virtue of the hydrodynamic hypothesis (|5.2p . Now it is possible to close the equations between den- 
sities and currents. Using again the hydrodynamic hypothesis with the fact that G is a. smooth deterministic 
operator, (|5.5p leads to, 

/ N Giijx + y,x)pa{x) + Qi2{x + y,x) 

Pa\X + y) — n r I \ 

y22[x + y,x) 

By differentiating this last relation w.r.t. y, altogether with l|5.6p and (|5.7p . we obtain the final deterministic 
expression for the current 



Ja{x) = D[ —2- + ^ aakPkPa 



k=2 



(5.8) 



which, combined with the continuity equation 



dpa ^ dJg 
dt dx 

leads again to a Burger's hydrodynamic equation. 



= 0, 



[The hetero-difFusion case] Here, the limit l|5.4p is a bit more tricky. In fact, the expansion of S involves 
correlations between currents and densities which already appear in the leading terms, and we expect an 
effective diffusion constant of the form 



n 



Da{p)=De^p(j2l3"''Pk 

k=2 



with 



ak def 



lim log 

N~too 



A, 



k=2 



k 



We pursue no further the study of this case, which presumably could be handled with block-estimates 
techniques (see [30]). 



5.2.2 DifFusion with reaction 

Here we treat the square-lattice model, a special case of (|2.2p .where reactions take place, in addition to 
diffusion. The procedure follows the lines of the preceding subsection. Using the mapping (|2.5p . the model 
is formulated in terms of two coupled exclusion processes, and the current equations corresponding to both 
species have the form 

Jr{^) = Kmr^+i ~ A-(z)ffrf+i, 
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with the rates given by l|2.6p . and we restrict the present analysis to the symmetric case (see relations (|2.4p ). 
Reversing for example the equation for Ja leads to the homographic relationship 

, _ AjWrf-jr'(i) 



which again can be iterated by means of a 2 x 2 matrix product, after defining u° and vf s.t. rf — wf/wf , 
Vie {I...N}. Define 



and 



7(iV) 



def T+W+T^W 



Then the proper scalings for large N are given by 



Letting now, 



\(N) 
lim = £>, 



lim . r, 



iim = 77. 



A^W _ /AIM 

Aa(») V K(^) 







G cannot be given explicitly, it is instead solution of the following combinatorial self-consistent equation 
The iteration now reads, 



H+p+l 



Vi+p+1 



G 



■ i+p i ^ 



Vi 



For the same reason as before, the limit Q oi G when N ^ go does satisfy 



x+y 



g{x + y,x) = G^{x + y,x) + / du g°{x + y,x + u)T,{x + u)Q{x + u,x), 



(5.10) 



with 



g°{y,x) = exp(j]cr I {2pb{u) - l)duj , 



by just applying the law of large numbers in the formal expansion of G with respect to E. We leave aside the 
question concerning existence and analytic properties of a solution of (|5.10p . We must again discriminate 
between two situations. 



[Case 7 = A] 



E(x) 



vi^Pb - 1) 



D 



_27/(2pb-l) 77(1 -2pb) 
which leads to the following differential system, 

^=77(2p,-lK--J„(a;)7;«, 

dv"- 

— = 27,(2pb - IK + ?7(1 - 2pbK, 
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5 Transient regime and fluctuations 



after making use of the law of large numbers and the hydrodynamic hypothesis. Combining these last two 
equations to express pjj — {u'^Va — v'^Ua)/v'^, leads to the relation 



[Case 7 ^ A] 

Like in the hetero- diffusion case of the last section, the effective diffusion constant Da{p) involves correlations 
between and Tf_^i and Ja{i/N) in the leading order term, and we expect a behavior of the form ^14j 

Da{pb) = Dexp^2pb{l - pb) log - 

as a result of a multipHcative process. This could be obtained through renormalization techniques appHed 
directly to equation l|5.9p . 



To conclude this section, we see that, for 7 = A, the differential system expressing, at steady state, the 
deterministic limit of the square lattice model with periodic boundary conditions finally reads, setting = 
2pa.b - 1, 

= 77(1 - vl)vb + Wa + 

t (5-11) 

— = -77(1 - I/j, jl^a + Wb + , 

where w is a possibly finite drift velocity and — ^p{va, Vb) and f''{i^a, i^b) are two constant currents in the 
translating frame. These currents have to be determined in a self-consistent manner, after fixing the average 
densities 7?a and Db and the periodic boundary conditions. For v — 0, the system (|5.1ip is Hamiltonian with 

= I ["a'^b -^l- ^b] + fbya - '-PaVb- (5.12) 

Indeed, it is easy to observe that IjS.lip can be rewritten as 

dva dH dvb dH 
dx dvb ' dx dva 

The degenerate fixed point i>a.b{x) — i'a,b is always a trivial solution and corresponds to the relations 

^ -qivl ~ l)vb, ipb ^ -nil ~ D^)Da- 



5.3 Microscopic currents 



5.3.1 Particle currents 



An important feature of our particle systems is that the number of particles is locally conserved. This 
property is refiected as iV ^ oo by a continuity equation, which relates local variations of particle density 
to inhomogeneous currents. In a discretized framework, conservation of particles is expressed according to 
the following 

Proposition 5.1. Let {Ji{t, e)} i — 1, . . . ,N be stochastic variables corresponding to the fl,uxes of particles of 
type /c e {1, . . . , Ti} between site i and i + 1, such that 

J^{t, e) i ^(Xf (t)X,'+i(t)Xi(i + e)X^^,{t + e)- X\{t)Xl,{t)X': {t + e)Xl,{t + e)) 

l^k 

with e > 0. By definition J^ {t,e) are ternary variables in {— i, 0, The following identity, equivalent to 
particle conservation, 

lim ^'^^ + ~ ^'^^^ + Jf^i(^, e) - (t, e) = a.s., (5.13) 
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holds for all i E {1,...,A^}, Vt G R+. In addition, letting ri^"\t) denote the sequence {X^{t)}, i = 
1, . . . ,N;k — 1, . . . , n}, then the variables {Ji(t, e)}, i — 1, . . . ,N;k ~ 1, . . . , n}, have a joint conditional 
Laplace transform given by 



,N def I 



C(0) 



k<l 



1) + XikX^X^^^ 



i=l 



where 4>k, k — 1, . . . ,n is a set of C°° bounded functions, and ipki — 4>k — 4>i- 



oie), 



(5.14) 



Proof. The points are mere consequences of the Markovian feature of the process and of its generator. In 
particular, (|5.13p results from the fact that, almost surely, at most one jump takes place in the time-interval 
e, when e — > 0, since all events are due to independent Poisson processes. In addition, on the time interval 
[t,t + e], the occurrence of a particle exchange between sites i and i + 1, corresponding to €J^{t,€) = 1 is 
only conditioned by the presence of a pair {k,l) at + 1), with a transition rate given by XkiXfX\j^^. 
Therefore 

i = l 

which, after a first order expansion with respect to e, leads to l|5.14p . 



5.3.2 An iterative numerical scheme 



Given a sample path r7'"'(i) at time t, we may generate a current sequence {Ji{t, e)} according to the local 
product form encountered earlier. In turn, once the set {Jf (t, e)} is known, the sequence r]{t + e) is almost 
surely determined, as e — > 0, by the identity (|5.13p . expressing conservation law of particles. We therefore 
have at hand an explicit stochastic numerical scheme to generate the sequence rjlt) step by step. 

Proposition 5.2. For any e > 0, G N, the iterative scheme given by 

ri' 

where Pe{ri\ri') is defined according to f5.13\) and ^5A^, converges when e ^ to the original probability 
measure Pt=ne (v) corresponds to the original process. 

Proof. There is only one thing to show: VT > 0, the probability Pe that 3t £ [0,T], such that two adjacent 
transitions occur within the same time-interval [t,t + e], tends to when e ^ 0. This is warranted by the 
fact that the total number of transitions for t < T is almost certainly finite. Indeed, we have 

NT 

^ kl ' e^O 

For the hydrodynamic limit the rates A^; scale like A^^ for large iV. Thus, it will be convenient to take a 
smg le limit e = e(Ar) ^ as iV oo, since the condition for the scheme to be meaningful writes 

iVe(Ar)(maxAfeO' = o(l), 

kl 

SO that we get a scaling of e(A^) = o{N~^) to meet our needs. This will allow us, in the sequel, to make use 
of the approximation 

N 

Y.a\(x^{t + 6) - x^it) Y.iJ'-i - J')') = 

i=l I 

for any set of bounded complex numbers {aj'}. ■ 



30 



5 Transient regime and fluctuations 



5.3.3 Central limit theorem for currents 



We are in position to exploit the conditional product form l|5.14p to perform a mapping, in the spirit of 
Lemma 4.1 of [14] , allowing to obtain a dynamical description of the system, in terms of some external free 
random process. To this end we assume, as a basic point, the hydrodynamic limit holds and we rest on the 
following lemma. 

Lemma 5.3. Suppose the existence of a set of density 

n,N . n,N 



E 



Pk, such that 



k=l 



for any given bounded complex function (j)^, and let 4> — sup {(j)k{x)). Then, 



k€{l..-7^} 

xe[o,i] 



E 



N . . ^ N . . . . , 



k<l 
i = l 



k<L 



From this we deduce the following identity, 

n,N 

,1 + 1 



C'(^) = exp(e^ A,,p'=(^)p'(^)(e^'^-(*) - 1) + W(-^)p'^(^)(e-^'^-(^) - l)) ) , 

k<l / 

(5.15) 

which leads to recover (in our specific context) a formulation of the general result of [3] concerning fluctuation 
laws of currents for diffusive systems. 



Keeping up to quadratic terms w.r.t. to functions (/)'s its argument, /ij"'((/)) reads, 

n,N . . „ n . . . ,2 

C(0)=exp(.^0.(-^):^'=(p(^)) + ^ ^ Mi^)Q^^iij)Hi^)+oi^)) (5.16) 



k=l 
i=l 



where are deterministic currents expressed, in terms of densities, by 



and Q is a n X n symmetric matrix 



OiklPkPl 



'h ^ Pi{l~ Pi). 

Q is of rank n — 1, because due to the exclusion constraint, currents are not independants, 

n 

^Jf(t,e) = 0, Vze{l,...,7V}. 

fc=i 

Let M the reduced matrix obtained from Q by deleting last row and last column. Its determinant is YYk=i Pk, 
so that it M invertible if none of the pk vanishes, with 
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after having taken into account the exclusion condition X]fc=i Pk — ^- Since every Hne k or column k sums 
to PkPn > 0, all the eigenvalues are strictly positive, and hence M{p) ows a real square-root matrix M^{p). 



Proposition 5.4. Let i 



1 denote a set of C°° hounded functions of the real variable x £ [0, 1], 



{<,fc = 1, 



1} a set of independent identically distributed Bernoulli random variables with parameters 



1/2, taking at time t values in { — 1/2, 1/2}. Then there exists a probability space, such that 



n,N 

k=i 
i=i 



1 

N 



E 

k=i 
i=i 



i^kn{^)[jHp{^))e + V2D'eJ2Mt{p{^)) ^^]+0(A^-'), a.s., (5.18) 



1=1 



The lines of arguments bare some features in common with the ones proposed in [14j (to study fluctuations 
at steady state). Recall, by law of large numbers, that correlations are negligible and do not affect the 
expression of the deterministic currents l|5.8p . This justifies the mapping IjS.lSp . On the other hand, the 

calculation of coefficents M-^- is done by comparing /ij"' in (|5.16p with 



E 



^ -1,N . 

exp(^ E ^-(]^)^^i(^)-O]=exp(^^0,(^)Q,(^)<^,(^) +,(,)), 



'^ = 1 kl 



because is symmetric and 

n—l . . n . . 

kl kl 

Setting, for A: = l,...,n— 1, 

[xN] 



I— 1 



the corresponding space time white noise processes 



W^(x,t) = lim — ^(x,t), 

W-»oo dx 



describe current fiuctuations in the continuous limit. 



5.4 Macroscopic fluctuations 

Two main quantities with be explored in this section: the Lagrangian and the large deviation functional. 



5.4.1 The Lagrangian 

The preceding section provides us with all coefficients required to achieve an informal derivation of the 
Lagrangian [3] describing the current fiuctuations. Given the empirical measure 

i=l 

and assuming the system admits a hydrodynamical description in terms of a density field pk{x,t), the 
statement in ||3j says that there is a large deviation principle for the stationary measure. In other words, the 
probability that the measure deviates from the hydrodynamic density profile pk is exponentially small 
and given by 

P{p<"'(t) ~ p{t),te [ti,t2]} ~ e-^^i'i.*^)(^\ 
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5 Transient regime and fluctuations 



where ^ 

%.t2](/5)= r L{p{t),dtm)dt. 
Jti 

Here the deviation from hydrodynamic solutions is due to current fluctuations. 

def X 

Writing = ^ the quantity V~^~9t — ^■^''{p)^ represents the fluctuations of the current . Reversing 
the relationship between current fluctuations and white noise process leads formally to 

-n-l 



^'^''^^^^^2UNT.^hJ'{^-'^+Mp)), ^-l,...,n-l. (5.19) 
' fc=i 

Then, replacing l|5.19p in the joint distribution of {W^{x^ t); x G [0, 1], A: = 1, . . . , n — 1} , we obtain 

c{m,dtm)dt = - / dxdt J2 {w\x, t)) 

k=l 



where e has been identified with dt and dx with Then, the symmetry of M 2 ^ the form l|5.17p of M 

and the exclusion constraint 



lead to the final compact form 



fc=0 



1 /-i -iv-'^ + Jkip) 



1 /■ 

•^0 J — 1 



k=l 



5.4.2 Hamilton-Jacobi equation and large deviation functional 

Here we proceed as in |3j. Let tt^, the conjugate variable of pk, 

, def dC{p, dtp) 

TTk[X,t) = — -r. 

ddtpk{x,t) 

The Hamiltonian is then given by 

f 1 ^ 



r " 

'^({Pfe,7i'fc}) / dx TTfc (a;, t)dtPk {x,t) - C. 
•^0 fc=i 

Algebraic manipulations lead to the expression 

7Y({pfe,7rfc}) = J dx dxTTkJkip) + Dpk(dxTTk^ 
Then the large deviation functional T, satisfying 

might be derived as in [3], from the following regular variational principle 

T{p) = inf /[_oc^o](/5), 
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where the minimum is taken over all trajectories p connecting the stationary deterministic equilibrium profiles 
Pk to p. This means that T and the action functional / must satisfy the related Hamilton-Jacobi equation 



In addition, one can check the relation 
where 



= 0. 



U= dx 



^akipk{x)pi{y)dy, 



•'^ k^i 

pi n 

S = - rfa; pfe log Pk, 
•^0 k=l 



a form already encountered in the reversible case, see equation p.iop . Indeed, when the process is reversible, 
U is translation invariant (i.e. independent of the initial integration point, here set to zero), and so 



dx- 



Jk 



' dpk{x) Dpk 
This approach could be used to analyse the non-reversible case. 



6 Concluding Remarks 

In this report we strove to put forward some techniques, and to extend methods to tackle the problem 
of mapping discrete model to continuous equations. Even in the context of a very specific model, namely 
stochastic distortions of discrete curves, some difficult questions remain. 

• The determination of the invariant measure in the general case, at the discrete level, which would 
generalize the totally asymmetric case [l7l [2^ . 

• The analysis of Hamilton-Jacobi equations to obtain a kind of continuous counterpart of the invariant 
measures, namely large deviation functionals. 

With regard to hydrodynamic limits, there is a puzzling issue, namely when particle-species diffuse at various 
speeds, in what we called the heterodiffusive case. For many one-dimensional models, it is well known that 
a single slow particle may considerably modify the macroscopic behavior of the system (see e.g. [25]) ■ For 
the time being, our approach is restricted to diffusive one-dimensional systems. Yet, other scalings (like 
Euler), as well as processes in higher dimension, are definitely worth being studied, in particular, it could 
be interesting to deal with more realistic exclusion processes, for instance those encountered in the field of 
traffic modelling. Besides, the analysis of irreversible invariant states in terms of cycles in a state-graph 
might well be extended to tackle ASEP on closed networks. 
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